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Abstract. We have considered the thermal equilibrium in pre-protostellar cores in the approximation where the 
dust temperature is independent of interactions with the gas and where the gas is heated both by collisions 
with dust grains and ionization by cosmic rays. We have then used these results to study the stability of cores 
in hydrostatic equilibrium in the limit where thermal pressure dominates over magnetic field and turbulence. 
We compare the density distribution derived in this manner with results obtained in the isothermal case. We 
find that for cores with characteristics similar to those observed in nearby molecular clouds, the gas and dust 
temperatures are coupled in the core interior with densities above ~ 3 x 10 4 cm -3 . As a consequence, one expects 
that the gas temperature like the dust temperature decreases towards the center of these objects. However, the 
regime where gas and dust temperatures are coupled coincides approximately with that in which CO and many 
other molecular species deplete onto dust grain surfaces. At larger radii and lower densities, the gas and dust 
temperatures decouple and the gas temperature tends to the value expected for cosmic ray heating alone. The 
density structure which one computes taking into account such deviations from isothermality are not greatly 
different from that expected for an isothermal Bonnor-Ebert sphere. It is impossible in the framework of these 
models to have a stable equilibrium core with mass above ~ 5 Mq and column density compatible with observed 
values (Ah > 2 x 10 22 cm -2 or Ay > 10 mag.). We conclude from this that observed high mass cores are either 
supported by magnetic field or turbulence or are already in a state of collapse. Lower mass cores on the other 
hand have stable states where thermal pressure alone provides support against gravitation and we conclude that 
the much studied object B68 may be in a state of stable equilibrium if the internal gas temperature is computed in 
self-consistent fashion. Finally we note that in molecular clouds such as Ophiuchus and Orion with high radiation 
fields and pressures, gas and dust temperatures are expected to be well coupled and hence in the absence of 
an internal heat source, one expects temperatures to decrease towards core centers and to be relatively high as 
compared to low pressure clouds like Taurus. 
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1. Introduction 

It has long been presumed that during the formation of 
a star, 
tostar' 



isons of the gravitational, thermal, magnetic and turbu- 
lent energies of such cores show that all these quantities 



there is an intermediate phase in which the "pro- 



is at least approximately in a state of hydrostatic 

equilibrium or magneto-hydrostatic equilibrium (see e.g. 



are equal to within the (considerable) uncertainties (Myers 
Goodman 1988a| ,fl; |Myers et al. 199l| ; |Crutcher 1999| ) 



Shu ct al. 1987 ) . While this idea in origin was merely based 



on plausibility arguments, it has received support from 
the discovery that in least some nearby molecular clouds, 
one finds embedded "cores" of higher density than the 
surroundings where the observed linewidths are thermally 
dominated. That is to say, while there may be local sub- 



These data suggest though they do not prove the existence 
of equilibrium structures which are an intermediate state 
in the evolution of a prestellar core. 

This picture has received more observational support 
with the advent of high quality maps of the millimeter 
continuum emission of the dust grains within such cores as 
well as the possibility to study their extinction in the near 



sonic turbulence, Lhei 



U IS 11U uv 



frtence for collapse unLu 



point mass and il appears LhaL Lheimal pressure is ca- 

pable of balancing gravity. On the other hand, compar- 



and mid infrared (lohnstone et al. 2000, 2001; Bacmann 



ct al. 2000| ; |Lada et al. 1994| ) . These have allowed a more 
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unbiased view to be obtained of the density distribution 
in prestellar cores. In particular, they have shown that 
the early molecular line data was highly biased because 



above a critical density of ~ 5 x 10 cm 3 (see Tafalla ct 
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al. 2002; |Caselli et al. 2002a| ,|b|; |Kramer et al. 1999| ) most 
molecular tracers including CO condense out onto dust 
grain surfaces. This depletion has the consequence that 
in molecular lines one sees mainly a lower density outer 
shell, whereas the dust emission (or absorption) offers a 
more unbiased view of the density distribution. 

The results from the millimeter continuum and in- 
frared absorption studies have been compared with a vari- 



ety of theoretical models of hydrostatic cores (Bacmann et 
al. 200f)|). One result of such studies has been evidence for 
a "flattening" in the density distribution for radii below a 
critical value of r cr ~ 2000-8000 AU implyi ng a steeper fall 
off in radius above r rr . Thus for example, Bacmann et al 



(2000|}]model L1544 with a roughly uniform density inside 
1900 AU but a rapid fall off outside this radius. Structures 
whose support against gravitational collapse is mainly due 
to the magnetic field are plausible both because the ob- 
served dust continuum maps show large departures from 
spherical symmetry and because of time scale arguments. 
Collapse on a free fall time scale would produce a larger 
star formation rate than that observed. 

There are however some cores where thermal pressure 
may dominate magnetic pressure and spherical symmetry 
may be a good assumption. Particularly worthy of note is 
B68 where Alves et al. (2001) have recently demonstrated 
that the density distribution derived from their NIR mea- 
surements is consistent with a purely thermally supported 
hydrostatic model. In fact, they find results consistent 



with the equilibrium structures discussed by Ebert (1955) 
and Bonnor (1956), which we will refer to in the following 
as Bonnor-Ebcrt spheres. To what extent B68 is excep- 
tional is presently unclear but we note that a large number 
of cores i n the Ophiuchus and Orion clouds ( Johnstone et 
al. 200[], 2001 ) appear to have characteristics compatible 
with Bonnor-Ebert spheres at temperatures of 15-30 K. 

This accord between theoretical expectation and ob- 
servation suggests that a study of the theoretical assump- 
tions may be worthwhile. One of these assumptions has 
been that of isothermality which is often based on the 
concept of "low optical depth for the cooling transitions" . 
In reality, the gas in prestellar cores is tho ught to be cooled 
mainly by optically thick CO transitions ( Goldsmith 2001 ) 
although, as mentioned above, the CO seems to disappear 
at high densities. Moreover, the typical density for which 
dust grains and gas become thermally coupled is roughly 
of the same order as that observed in prestellar cores (see, 
e.g. Kriigel & Walmsley 1984). Thus, we decided that a 
new look at the gas temperature distribution to be ex- 
pected in such cores seemed warranted. 

In this paper we replace the assumption of an isother- 
mal gas with the more realistic condition of thermal bal- 
ance in the gas, and we evaluate the consequences of an 
external heating source (the interstellar radiation field) 
on the structure of the cloud (density and temperature 
profiles) and its stability properties. Recent st udies of the 



dust temperature distribu tion in such objects ( Evans et al 
2001; ^ucconi et al. 2001 ) have shown that the dust tem- 



It seems reasonable to ask how the gas temperature will 
react in such circumstances. One might also ask whether 
the density distribution in hydrostatic equilibrium will de- 
part appreciably from that expected under the isothermal 
assumption. Will the expected density contrast differ for 
example from that expected for a Bonnor-Ebert sphere 
when one calculates the gas temperature in self-consistent 
fashion? This article represents an attempt to answer such 
questions. 

The outline of this paper is as follows. In Sect. 2 we 
give a brief introduction to the theory of structures in hy- 
drostatic equilibrium including the results for an isother- 
mal equilibrium Bonnor-Ebert sphere. In Sect. 3 we dis- 
cuss the input to our calculations and the simplifications 
which we have made. In Sect. 4, we present our results 
for the gas temperature distribution in two model cores 
whose density distribution has been assumed similar to 
that observed in L1544 and B68. In Sect. 5 we present 
our results for non-isothermal hydrostatic equilibria for a 
variety of conditions. Here we show among other things 
that the equilibria obtained depend relatively sensitively 
on the external radiation field. In Sect. 6 we discuss the ob- 
servational consequences of our results, and compare the 
properties of our model clouds to isentropic polytropes. In 
Sect. 7 we summarize our conclusions. 



2. Theoretical background 

An isothermal gas sphere embedded in an external 
medium of given pressure p cx t has a critical Bonnor-Ebert 
mass above which no state of equilibrium can exist ( Ebert 
1955| ; [Bonnor 1956] ; |McCrea 1957b , 



M BE - 1-182 



xj G 3 p cxt 



2.6 



/ T 



Pcxt 



\10KJ V 2xl 4 Kcm- 



-1/2 



(1) 



where a = (kT g / pmu) 1 ^ 2 is the sound speed in the gas, 
T g the gas temperature and p the mean molecular weight 
(see Fig. |). 

Below this critical value, the equation of hydrostatic 
equilibrium admits single or multiple solutions, charac- 
terized by different degrees of density concentration. It 
has been shown by Bonnor (1956| ) and Ebert (1957) that 
equilibrium configurations with center-to-boundary den- 
sity contrast p c /Pb < 13.98 are stable, whereas equilib- 
rium solutions with p c /pb > 13.98 are unstable to radial 
collapse. Stable isothermal equilibria can also be charac- 
terized by the condition 



R 



£ max = —y/AnGpe < 6.451, 
a 

where R is the cloud's radius. 



(2) 



perature typically falls by a factor of 2 from edge to center. 



Recently, |Lombardi fc Bcrtin (2001| ) have re-examined 
the classical Bonnor-Ebert problem relaxing the assump- 
tion of spherical shape of the cloud. The result is that 
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for non-spherical isothermal clouds, the critical Bonnor- 
Ebert mass is larger than the value given by Eq. (0); a 
sphercal shape is more prone to gravitational instability 
than clouds of other shapes. The density contrast p c /pb 
for marginally stable non-spherical clouds is different from 
the Bonnor-Ebert value 13.98 for spherical clouds, but, 
remarkably, the condition for stability, when expressed in 
terms of the ratio of the average density p to the bound- 
ary density pb, i-e. p/pb < 2.465, is independent of cloud 
shape. 

To express the condition for stability (M < Mbe and 
Pel ' Pb < 13.98) in terms of observable quantities one can 
eliminate the external pressure writing p oxt = a 2 p^ > 
13.98 _1 a 2 p c , obtaining, for a given cloud's mass M, 



p c < 19.53- 



(3) 



G 3 M 2 ' 

that, for a cloud of molecular hydrogen with mean molec- 
ular weight p, = 2.33 becomes 

" ' L<i ' '"' 1 1 (jT^ ''• (4) 



10 K 



Thus, stable isothermal equilibria are only possible pro- 
vided the central density is relatively low (in the ab- 
sence of magnetic fields): this causes a conflict with ob- 
servations if one attempts to model observed molecular 
cores with such equilibrium configurations. For example, 
if dense cores with M > 4 M@ and T g = 10 K are modeled 
as stable Bonnor-Ebert spheres, their density according to 
Eq. (U) will fall below the critical value of ~ 10 4 cm" 3 re- 
qired to excite the molecular lines observed in these re- 
gions. Moreover, the central densities inferred in many 
pre-stellar cores based on the observed mm-submm dust 
emission (e.g. Ward-Thompson ct al. 1999) is often above 



10 cm~ J and hence for gas temperatures of ~ 10 K, only 
cores of mass below 1 Mq can be in stable equilibrium. 

Many of the observed objects however have higher 
masses (see the discussion of Bacmann et al. (2000) and 



also column densities of up to lO^cm -2 in molecular hy- 
drogen or 100 visual magnitudes of extinction. One can 
express the stability condition for Bonnor-Ebert spheres 
as a condition on the cloud's total (from edge to edge) ex- 
tinction through the center Ay,tot = iV(H 2 )/10 21 cm~ 2 , 
obtaining 



.4 



V.tot 



^29 



IV 



10 K 



M 



mag. 



(5) 



This is clearly much less than often observed and provides 
a strong argument for the support of many observed struc- 
tures being due to forces other than pure thermal pressure 



(though there are also cases of thermal support, e.g. Alves 
et al. 2 j00l] ; [Johnstone et al. 2000| , |2001| ). In the following 
we consider "observed cores" as having a minimum A\ of 
10 mag., and thus exclude processes such as UV photo- 
electric effect from consideration. We do not moreover in 
this work compute models including magnetic fields and 
turbulence, but briefly mention in the next Sections how 
these agent of support modify the picture outlined above. 



In summary, non-magnetic stable isothermal cores rep- 
resent only moderate enhancements (a factor ~ 2.5 in the 
average density) of the ambient gas density, quite inde- 
pendent on cloud shape and the turbulent field. This lim- 
ited density range represents a serious limitation when 
modeling molecular cloud cores as non-magnetized hydro- 



static equilibria. For example Boland & de Jong (1984) 
were able to reproduce column density, brightness tem- 
peratures and abundances of several chemical species mea- 
sured in dark clouds like L134, L183 and TMC-1 only with 
models characte rized by density contrast p c /pb — 2000. 
Hascgawa (198§| ) however verified that these models were 
well into the unstable regime. The same cr iticism was 
made by Chieze fc Pincau des Forets (1987) to the hy - 



drostatic models computed by Falgarone fc Puget (1985 ) 
|Hasegawa (1988 ) concluded that stable spherical, hydro- 



static models able to satisfy all observational constraints 
on extinction, external pressure, mass, and chemical com- 
position of molecular cloud cores, could no be constructed. 
In this respect, it is interesting to recall that when molec- 
ular cloud cores are fitted by Bonnor-Ebert spheres, the 
best fit is obtained for models that arc in the unstable 
regime: £ max = 6.9 for B68 flAlves et al. 200l|), £ max = 12.5 



for B335 (Harvey et al. 2001, although B335 is not a star- 
less core). 

2.1. Magnetic fields 

The additional support provided by a large-scale magnetic 
fields can be evaluated by a virial-theorem analysis or with 



the help of detailed calculations (see e.g. McKee et al 



1987, for a review). The maximum value of a cloud mass 



for the existence of stable equilibria is approximately equal 
to 



M, 



1 - 



0.13 



-i -3/2 



Mi 



(6) 



where A is the magnetic mass-to-flux ratio of the cloud in 
units of G -1 / 2 ( Mouschovias fc Spitzer 1976j ; [Tomisaka ct 



|al. 1988 ). For magnetically subcritical clouds (A<0.13), no 
amount of external compression can induce the collapse of 
the cloud (as long as the magnetic field remains frozen in 
the matter). Magnetically supercritical clouds (A > 0.13) 
can exist in stable equilibrium for masses larger than the 
Bonnor-Ebert critical mass, and can reach values of p c /pb 
much higher than 13 .98 (up to ~ 100, see Fig. 4 and 6 of 
[Tomisaka ct al. 1988|) . 



2.2. Non-uniform turbulence and temperature 
gradients 

If the turbulent motions of the gas are characterized by a 
uniform (three-dimensional) mean square velocity (v 2 urb ), 
the effects of a turbulent pressure, in addition to the ther- 
mal pressure, can be trivially incorporated in the clas- 
sic Bonnor-Ebert picture summarized above. The condi- 
tion for the existence of stable equilibria is the same ex- 
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pressed by Eq. dll) with the sound speed a 2 replaced by 3.1. Dust and gas temperature 



an "effective" sound speed a 2 s 



a 2 + 



v turb)A allowing 



the existence of stable equilibria for values of the central 
density (or peak extinction) larger than for clouds sup- 
ported by thermal pressure alone. Notice, however, that 
(as in the purely thermal case) marginally stable clouds 
are characterized by a density contrast p c /pb = 13.98 or 
p/Pb = 2.465, where p is the average density, irrespectively 
on the cloud's effective temperature. 

One should remember also that the observed line 
widths (which in some cases are close to being thermal) 
place limits on possible values of v tur b- In dense cores, 
as well as in giant molecular clouds, turbulent linewidths 
are often observed to increase at larger scales (see discus- 
. Barranco & Goodman 1998|; |Goodman ct al. 1998) ) , 



Thus we determine the gas temperature T g by solving the 
equation of thermal balance of the gas in the presence of 
dust grains heated by the external radiation field, 



T rr - A e - A ffri = 0, 



(7) 



where r cr is the cosmic-ray heating rate, A g the gas cool- 
ing rate by molecular and atomic transi tions, and A E d the 
gas-dust energy transfer rate. We follow poldsmith (200l| ) 
and adopt a cosmic-ray heating rate 



r, r = lo- 



3 x 10" 



n(H 2 
cm" 



erg cm 



-3 _-i 



(8) 



sion m 



suggesting that these objects can be modeled as negative 
inde x polytropes (|dc Jong ct al. 1980 ; Dickman & Clemens 



1983; |laloncy 1988| ). This is easily understood compar 



ing the polytropic equation of state P cx p 1+1 / n with the 
expression for the thermal pressure of a gas, P cx pT, 
implying T cx p 1 /™. Since the density is always a decreas- 
ing function of radius, negative index polytropes represent 
equilibrium configurations in which the effective temper- 
ature increases outward. However, from the point of view 
of the stability properties, our models show a behavior 
similar to polytropes with n>l, We will return to this 
point in Sect. 6.2. 

The logatropic equation of state analyzed by 



The gas cooling rate A g for typical conditions of molecu- 
lar cloud cores has been computed by Goldsmith (200l[ ) 
in the LVG approximation (assuming a velocity gradi- 
ent of 1 km s _1 pc _1 ) for various degrees of depletion 
of molecules from the gas phase onto grain surfaces. We 
adopt the parametrization of Goldsmith (2001) 



A D 



10 K 



erg cm 3 s , 



(9) 



where a and [3 are parameters that depend on the H 2 
density and the molecular depletion factors. When not 
explicitly noted, we assume that molecular abundances in 
the gas have their standard (undepleted) values. As for the 
the gas- dust energy transf er rate we adopt the exp ression 



McLau ghlin fe Pudritg (1 -906 ) to incorporate the contribu 



tion or on non thermal motions to the aupport of molecu — 

lar clouds (and cloud cores) makes possible the existence of 
stable equilibria for large density constrasts (p c /pb — 100) 
but the resulting configurations are much more extended 
than Bonnor-Ebert spheres and are characterized by a 
maximum ratio p/pb — 3/2, even smaller than for crit- 
ical Bonnor-Ebert sph eres. Indeed, the obse rved column 
density distributions flBacmann ct al. 2000| ) of cores do 
not fit well to logatropes. 



3. Model assumptions 



Hollcnbach (1983| ) 



given by Goldsmith (2001 ) based on calculations of 



Burke 



A P 



2 x 10 



-33 



Vcm 13 / 



2 (T v -T* 



K 



10 K 



1/2 



erg cm 3 s . 



(10) 



To determine the dust temperature we ignore the 
gas-dust coupling and we solve the equation of thermal 
balance of the dust 



- A d = 0, 



(11) 



where r cxt is the dust heating rate by the external radi- 
ation field and Ad the dust cooling rate by radiation. In 



Followi ng the approach outlined by — Goldsmith (2001 ) 



practice, we use the analytical solution of Zucconi et al 



we co mpute the gas thermal equilibrium assuming that 



(200l[ ) to compute Td as function of extinction from the 



the dust temperature is unaffected by collisions with gas 
particles (see Evans et al. 2001, for estimates of the va- 
lidity of this). Moreover, we have used the analytical fit 



cloud's boundary. 
As in 



Zucconi et al. (2001), we have adopted as a stan- 



of Zuc qoni et al. (20 61 ) to dotormino dust temperature as 



dard the interstellar radiation field (ISRF) in the solar 
neighborhood given by |Mathis et al. (1983j ) and |Black| 



a func iion of radius in our spherically symmetric model 

cores. The justification for this is simplicity. We are inter- 
ested here in understanding qualitatively to what extent 
departures from isothermality may affect the properties 
of pre-stellar cores. High precision in determining the gas 
cooling rates and the dust temperature distribution is un- 
likely to fundamentally affect our results. 



(1994). It consists of four components: the V-NIR com- 
ponent, peaking at A ~ 1 pm, the MIR and FIR compo- 
nents, peaking at A ~ 100 and 140 pm respectively, and 
the cosmic background radiation, peaking at A ~ 1 mm. 
We allow for variations of the ISRF, assuming, for simplic- 
ity, that the spectral shapes of these components are the 
same as the local ISRF, and we scale the intensity of the 
V-NIR, MIR, and FIR components by a factor Go (with 
Go = 1 for the local ISRF). 
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The existence of a cloud radius leads to the necessity of 
an external pressure to keep the system in equilibrium. We 
embed our model cloud cores in a spherical envelope that 
represents the ambient molecular cloud. We do not need 
to specify the physical properties of this envelope; we just 
assume that it provides the needed external pressure and 
a shielding from the ISRF directly incident on the cloud 
cores. We set the latter at the value A™ v = 1 (McKee 
1999|)7|sufhcient to absorb the UV and most of the V- 



NIR radiation. We define Ay jC to be the centre-to-edge 
extinction of our model cloud cores. The range of masses 
and densities of interest here corresponds to the range 
5 < v4-v,c is 100, where the central dust temperature is 
determined mostly by optically thick absorption of V-NIR 
radiation and optically thin absorbtion of MIR and FIR 
radiation. 

The fundamental parameters of our models are the ex- 
ternal pressure p C xt arid the intensity of the external ra- 
diation field Go- In this paper, we consider independent 
variations of these two quantities, and we analyse sepa- 
rately the dependence of the critical mass for increasing 
p cxt and fixed Go (Sect. 5, Fig. ^|), and the effect of varying 
Go for a given value of p cxt (Sect. 6, Fig. ||). 

3.2. Effects of molecular depletion 

We include in our models the effect of molecular depletion 
onto dust grains, parametrized for simplicity by the single 
parameter /<j (depletion factor), defined as the ratio of the 
"standard" abundance of CO isotopes and CS to the ac- 
tual abundances. The corresponding abundances of other 
molecular species are given in Table 3 of Goldsmith (200 lD . 
Gas-phase cooling rates for several values of fa in the 
range 1-100 were computed and tabulated by Goldsmitl] 



We 



|(2001| ) for H 2 density from 10 3 cm" 3 to 10 6 cm 
assume the dependence of /d on density given by Fafalla 
|et al. (2002| ), 



/d 



exp 



where the critical density for CO depletion is taken to be 
«dc P = 5.5 x 10 4 cm -3 . 

Although we are mainly interested in applications to 
molecular cloud cores of moderate to high extinction, 
where molecular depletion is effective only in a region of 
limited spatial extent, one should keep in mind that large 
fraction of the cloud's mass is contained in the outer parts 
(roughly M(r) oc r), and therefore the overall stability of 
the cloud may be modified by the thermal properties of 
these external layers. 

In our models we find in general that the effects of de- 
pletion on the temperature of the gas are significant only 
for densities n < 10 5 cm" 3 . At higher density, gas-dust 
coupling overwhelms the reduction of the cooling rate of 
the gas, and the gas temperature becomes insensitive to 
depletion. This is in qualitative agreement with the re- 
sults of Goldsmith (2001 ). We conclude with GoldsmitT] 



(2001) that molecular depletion in the central, dense re- 
gions of molecular cloud cores (n> 10 45 -10 5 cm" 3 ) is not 
associated with an increase in the gas temperature. 

In the extreme case where CO (and other molecules 
containing C,N,0,S, etc.) depletes out onto grains causing 
Eq. (0) to reduce to T cr = A g d. Then using Eq. (||) and 
@, one finds 



< 



3 x 10" 



— )( 



. 10 4 cm" 3 



10 K 



-1.5 



(13) 



This shows that in the extreme case where all heavy 
species deplete out, coupling between gas and dust tem- 
peratures will be assured for densities larger than 5 x 
10 4 cm" 3 if the cosmic ray ionization rate is standard or 
lower. Since it is essentially at densities above this value 
that CO is observed to deplete out, we conclude that a 
rough equality between gas and dust temperatures in the 
depleted region is a reasonable assumption. In this con- 
text, it is interesting to note that Bergin et al. (2002 ) sug- 
gest that even such relatively volatile species as molecular 
nitrogen have depleted out in the central core of B68. 



4. Temperature distributions for L1544 and B68 

Before considering what one expects for theoretical equi- 
librium cores, it is useful to consider what one predicts for 
the temperature variation within an object whose density 
structure is similar to those observed. Specifically, we will 
consider the case of LI 544 flCaselli et al. 2002a|,|b|) and B68 



Alves et al. (2001), two well studied prestellar cores. 

To model L1544 we assume Go = 1 and the den- 
sity structure of Tafalla et al. (20*02] ), which, in turn, is 
based on the 1.3mm dust emission and absorption maps 
of [Bacmann et al. (200C| ). Thus, we assume the H 2 density 



distribution n(r) to be given by 



(12) n (r) = 



1 + (r/r ) 



2.5 ' 



(14) 



where the central H2 density is n c = 1.4 X 10 6 cm" 3 and 
tq = 0.014 pc. It is worth noting here that we have on pur- 
pose taken here a somewhat extreme case with large den- 
sity contrast relative to the surroundings and large central 
column density (6.0 x 10 22 cm" 2 corresponding to 60 mag. 



of extinction, from Bacmann et al. 2000). One should also 
be aware that the aspect ratio of the core based on the dust 
emission contours is ~ 0.6 and hence spherical symmetry 
is a crude assumption. Nevertheless, the case is of interest 
in that it is a model which approximates the actually ob- 
served density distribution and in which, one expects that 
the central temperature will be relatively low (due to the 
large extinction). 

For B68 we assume the density structure of a Bonnor- 



Ebcrt sphere with the physical parameters given by Alvc; 
|et al. (200 l"l ). The central density n c = 2.4 x 10 5 cm" 3 
is a factor ~ 6 lower than in the case of L 1554, and 
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consequently the amount of CO depletion at the center 
of B68 is about two orders of magnitude lower than at 
the center of L1544. The intensity of the IRSF incident 
on B68 can be estimated as follows. We first estimate the 
external UV radiation field based on the 90 /jm data of 
Ward- Thompson et al. (2002) and the conversion between 



FIR intensity and incident field found by Boulanger & 
Peraul; (1988). Another approach is to use the 7fim in- 
tensity from Bacmann et al. (2000) ) and the results of 
Boulanger et al. (1996) to infer Go- Combining these, we 
conclude that the radiation field incident on B68 is ~ 2.5 
times larger that the standard ISRF. 

In Fig. [|, we show the inferred gas temperature de- 
pendence for L1544 and B68, computed with cosmic ray 
ionization rate £ = 1.3 x 10~ 17 s -1 and accomodation co- 
efficient = 0.3. The assumed CO depletion factor (cut off 
at a value of 100) is shown for comparison. In the inner 
part of both cores (densities above 3 x 10 4 cm -3 ), the gas 
temperature is coupled to that of the dust and thus de- 
creases gradually towards the core center. This coupling 
is more effective in the case of LI 544, which has a higher 
central density. The gas is slightly hotter than the dust 
however due to cosmic ray heating and this difference in- 
creases as the density decreases. The CO depletion ap- 
pears to have only small effects on the gas temperature 
profile because while at high density, gas-dust coupling 
dominates, at densities below 10 5 cm -3 , the depletion is 
not very large (and the optically thick nature of CO cool- 
ing causes the dependence of CO cooling on abundance to 
be minor). We draw the conclusion from this study that 
in the depleted core nucleus (densities above 10 5 cm -3 ) 
where many molecular species may be unobservable be- 
cause they are in solid form, the gas temperature may be 
expected to be close to that of the dust. In outer regions 
where the density is below 10 5 cm -3 and where shielding 
from the external radiation field is much less, the gas and 
dust temperatures become uncoupled and one may find 
gas temperatures lower than that of the dust. 

It is worth pointing out that the scenario depicted in 
Fig, p] yields gas tem peratures in conflict with observation. 
Tafalla et al. (2002] ) find based on their NH3 observations 
a temperature of 8.7 K throughout the L1544 core with no 
evidence for density gradients (40" HPB W equivalent to 
0.03 parsec), whereas Bourkc et al. (1995) give T g — 16 K 
for B68, also derived from NH3 data. These are in rough 
agreement with the values in Fig. [l] though the LI 544 
data do not show evidence for increasing temperature with 
radius as predicted. 

Another conclusion that one can draw from Fig. [I] is 
that the gas temperature and hence pressure in the high 
density core nucleus depends somewhat on the external 
radiation field even for cores of extremely high visual ex- 
tinction. Thus one can expect cores in Taurus for exam- 
ple (where the incident radiation field appears to be aver- 
age) to have different characteristics than Ophiuchus cores 
where the radiation field appears in general to be an or- 
der of magnitude higher. Of course, there are other effects 
such as the external pressure which may differ in the two 




J 
20 




0.02 0.04 0.06 

radius (pc) 

Fig. 1. Upper panel: inferred gas (long-dashed curve) and 
dust temperature dependence (short-dashed curve) for the 
model of L1544 discussed in the text (with Go = 1). The 
assumed density distribution from Eq. ( |l4[ ) is shown for 
comparison (solid curve). The assumed CO depletion fac- 
tor /d (multiplied by 10 4 , dotted curve) is also shown. 
Lower panel: same as above, for the model of B68 (with 
Go = 2.5). One sees that for n > 3 x 10 4 cm -3 the gas 
temperature follows the dust temperature and decreases 
towards the core center. Outside this rigion, gas and dust 
are decoupled and the gas heating is determined solely by 
cosmic ray ionization. 

cases also but we stress that the radiation field alone may 
cause core temperatures and hence accretion rates to be 
higher in Ophiuchus than in Taurus. 

5. Self consistent hydrostatic equilibrium results 

The equation of hydrostatic equilibrium for a spherically 
symmetric cloud is 



dp + dV _ Q 
dr ^ dr 



(15) 



where p is the gas pressure and V the gravitational po- 
tential. The latter satisfies Poisson's equation, which in 
spherical symmetry reads 



1 d 



,dV\ 



AnGp. 



(16) 



r 2 dr \ dr J 

The gas pressure p, neglecting the contribution of turbu 
lent pressure, is 

9 



P 



-kT« 



(17) 



where T g is the gas temperature, computed as described 
in Sect. 3, and /1 = 2.33 the mean molecular weight. 
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The behavior of our model clouds with respect to an 
external compression is qualitatively similar to that of 
isothermal spheres. Fig. [2] shows the locus in the radius- 
external pressure diagram of cloud models with M = 
1 M©, for different values of the ISRF intensity Go. For 
large radii the cloud follows closely the T g = 10 K Bonnor- 
Ebert isotherm, independent of the radiation field (the ef- 
fect of gas-dust coupling is negligible for these diffuse con- 
figurations). At higher values of the external pressure the 
cloud radius is smaller and the central (and average) den- 
sity becomes larger. As a result, the central gas tempera- 
ture increases and the sequence of equilibria departs from 
the Bonnor-Ebert isotherm. Eventually, a critical value 
of the pressure is reached above which no equilibria exist, 
and the curve turns down. Past this critical point, equilib- 
ria are characterized by an increasing central density and 
density contrast, and are unstable (see following Section). 
For clarity, in Fig. ^ we truncate the curves of equilibrium 
slightly after the critical point. 

Fig. H shows the density of H2, the temperature of the 
gas and the temperature of the dust as function of radius 
for the marginally stable clouds of mass M — 1 Mq with 
G = 1 (central H 2 density n c = 2.7 x 10 5 cm" 3 ). The 
density profile is very close to that of an isothermal sphere 
with the same values of the central density, shown for com- 
parison in Fig. U for temperatures equal to the central and 
boundary temperature of the model. The gas temperature 
decreases from the outer regions to the center of the cloud. 
For marginally stable clouds with M — 1 Af Q the central 
density is sufficient to ensure a modest gas-dust coupling, 
but the gas temperature is mainly determined by the bal- 
ance of cosmic-ray heating and gas cooling. 

Fig. [| shows the H2 density, the gas and dust tem- 
peratures and the extinction profile of the marginally 
stable cloud in thermal equilibrium with an ISRF with 
Go = 100 and bounded by an external pressure p C xt = 
1 x 10 6 K cm~ 3 . These roughly corresp ond to the struc- 
tures observed by Johnstone et al. (2001 ) in the pOph and 
Orion clouds. For comparison, the figure also shows the 
density profiles of marginally stable Bonnor-Ebert spheres 
with the same central density and uniform temperature 
equal to the gas temperature at the center and at the 
boundary of the cloud. 

In Tables |] and || we list the values of central H2 den- 
sity, gas and dust temperature, density contrast, and op- 
tical extinction at the center for marginally stable clouds 
with M = 1 and 5 M Q and various scaling factors for the 
ISRF. In these models we have varied p ex t until we find 
the marginally stable model. 

A comparison of the results for the two different masses 
is instructive. Small mass clouds (M ~ 1 Mq) exposed 
to the average ISRF behave like isothermal spheres with 
T g ~ 10 K, and are characterized by a maximum density 
contrast and critical central density very close to those 
of the corresponding Bonnor-Ebert values. Increasing the 
ISRF (Go = 10) has the effect of raising the dust tempera- 
ture, and of increasing slightly the gas temperature in the 
central regions. This increase in the central gas tempera- 



ture causes a higher central density in the marginally sta- 
ble core and this, in turn, increases the efficiency of dust- 
gas coupling. A further increase of the ISRF (Go = 100) 
establishes the approximate equality of the gas and dust 
temperatures at the center. Thus, the fact that small 
mass clouds are dynamically stable for larger values of 
the central density than their more massive counterparts 
(see Eq. |^), has the consequence that they are "stabi- 
lized" by an external ISRF simply because they are able 
to become warmer if the dust does. Their gas tempera- 
ture profiles, however, remain quite uniform, and therefore 
the maximum density contrast cannot deviate significantly 
from the Bonnor-Ebert value p c / Pb — 13.98, as shown in 
Table @. 

A different behaviour characterizes more massive 
clouds (M ~ 5 Mq), as they become dynamically un- 
stable before an efficient gas-dust coupling can be estab- 
lished, for any value of the intensity of the ISRF (see 
Table ||). Their central extinction is low (Ay, c ~ 3-5), and 
hence they do not correspond to "real observed cores". 
Nevertheless, it is interesting that they can develop sig- 
nificant inward temperature gradients, because molecular 
cooling is less efficient than cosmic-ray heating in the cen- 
tral regions, characterized by densities n ~ 10 4 -10 5 cm~ 3 
(where A g oc n - 5 ) than in the outer parts where n ~ 10 3 
(and A g oc n, as cosmic-ray heating). An inward tempera- 
ture gradient stabilizes the cloud (see Sect. 6), and allows 
large deviations in the maximum density contrast p c /ph, 
up to 30-40, as shown in Table |[ Notice that any molec- 
ular depletion effect, ignored in these calculation, would 
make the central temperature even higher, and thus in- 
crease the slope of the temperature gradient. 

In Fig. [5] we summarize our results for the maximum 
density contrast p c /pb allowed for clouds in the mass 
range M — 1-10 Mq. Again p ext has been varied until 
the marginally stable model is found. Deviations of p c /pb 
from the Bonnor-Ebert value 13.98 indicate deviations 
from a condition of uniform gas temperature in the cloud. 
As anticipated in the above discussion, the largest devia- 
tions from isothermality are obtained for masses M ~ 4- 
6 Mq. Marginally stable clouds of small mass (M<2 Mq) 
are characterized by relatively large values of the cen- 
tral density, and by a sufficiently good coupling of gas 
and dust: they are uniformly warmer or colder accord- 
ing to the intensity of the external ISRF, but because of 
this coupling they do not develop significant temperature 
gradients, and therefore behave essentially as isothermal 
spheres. Marginally stable high-mass clouds (M > 8 Mq) 
also are similar to isothermal spheres, but for a different 
reason: they are characterized by relatively low values of 
the central density (n c ~ 10 3 cm -3 ) and also have low 
Av, c (unlike observed cores), which implies an uniform 
value of the gas temperature since both cosmic-ray heating 
and radiative gas cooling depend approximately linearly 
on density (and gas-dust coupling plays no role at these 
densities). In the intermediate mass range, both because 
of the different dependence of heating and cooling rates 
from density and partial gas-dust coupling, the deviations 
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Fig. 2. Sequences of model clouds with M = 1 Mq and 
Ay™ — 1 in the radius-pressure diagram, for various scal- 
ing factors Go of the ISRF (solid lines). For compari- 
son, the dotted line shows a Bonnor-Ebcrt isotherm with 
T g = 10 K for the same cloud mass. 

Table 1. Properties of marginally stable clouds 

WITH M = 1 Mq 



Fig. 3. Density of H2 (solid line), temperature of the 
dust (short-dashed line), temperature of the gas (long- 
dashed line), and visual extinction (dash- dotted line) as 
function of radius for the marginally stable clouds with 
M = 1 M Q and Go = 1. The dotted lines show for com- 
parison the density profile of an isothermal sphere with 
the same values of n c and T = T g c or Tg.b- The visual 
extinction at the cloud's center is Ay.c = 14. 
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Table 2. Properties of marginally stable clouds 
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from uniform gas temperature are larger, and the maxi- 
mum density contrast p c /Pb na s a peak for M ~ 4-5 Mq. 
Fig. ||also shows, for comparison, the values of the density 
contrast for marginally stable polytropic spheres of index 
n, discussed in the next Section. 



10" 
10 




0.01 
radius (pc) 



Fig. 4. Same as in Fig. |, but for a marginally stable cloud 
bounded by an external pressure p cxt = 1 X 10 6 K cm -3 
and with G = 100. The mass of the cloud is M = 1.6 Mq 
and the extinction from the center to the edge of the cloud 
is A v = 29. 
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mal spheres. According to the static method, continuous 
(dashed) lines represent stable (unstable) equilibria. The 
dotted lines indicates the boundaries of the stability do- 
main for isothermal spheres with arbitrary temperatures 
(only the case T gas = 10 K is shown) and non-isothermal 
spheres in thermal equilibrium with the external ISRF. 
The figure shows that for a given external pressure the 
condition of thermal equilibrium allows the existence of 
stable cloud configurations for values of mass and central 
density not permitted under the condition of uniform gas 
temperature. 

Although not very large, this effect may resolve in part 
the problem posed by objects like B68 which are well fitted 
by slightly unstable Bonnor-Ebert spheres (magnetic fields 
are of course another possibility). For example, for the 
three marginally stable non-isothermal clouds of Fig. |^ 
the central-to-boundary density ratio are (p c / Pb) = 19; 
25 and 39, for Go = 1, 10 and 100, compared with the 
corresponding ra tio 13.98 for isothe rmal spheres, and 16.5 
inferred for B68 flAlves et al. 200l|). 



Fig. 5. Value of the central-to-boundary density ratio 
(Pc/Pb)cr as function of the cloud's mass for different scal- 
ing factors Go of the ISRF (dots and solid lines). The 
dashed lines show the corresponding values of (p c / pb)cr = 
13.98 for marginally stable isothermal and polytropic 
spheres of index n. 



6. The effect of the external radiation field 

The presence of an external ISRF affects the stability of 
clouds bounded by an external pressure. In the first place, 
it makes the dust (and therefore the gas) warmer on av- 
erage, thus increasing the critical mass for gravitational 
instability (proportional to Tg ); in the second place, it en- 
forces a temperature gradient in the cloud which modifies 
the value of the classic Bonnor-Ebert critical mass. 

The latter point is best investigated by the so-called 
static method originally devised by Zeldovich (1963) for 
white dwarfs and neutron stars (see Tassoul 1978). This 
method allows one to perform a stability analysis for ax- 
isymmetric motions on the basis of the properties of the 
equilibrium models only. The static criterion asserts that 
a turnover in the M(p c ) curves at fixed P ex ti occurring at 
say p c — Pc r (Pcxt) ; marks the onset of dynamical insta- 

= 



bility to radial motions. The condition (dM /dp c ) Pc 
is both a sufficient and necessary condition for a stability 
transition: equilibria with p c < /?J: r (pext) are stable, while 
equilibria with p c > p"{p cx t) are unstable to at least one 
normal mode. 

In Fig. H we show how the intensity of the external 
ISRF may affect the stability of a cloud bounded by a 
given external pressure. For this part icular examp le, we 
have chosen p cxt = 2 x 10 4 K cm" 3 ( |McKee 1999j) . The 
three curves labelled with Go = 1, 10 and 100 show 
the mass versus central H 2 density for externally heated 
clouds, compared with the similar quantities for isother- 



6.1. Comparison with polytropic models 

From the point of view of the stability properties and the 
maximum density contrast, the models shown in this pa- 
per have a behavior similar to polytropes with b>1 (see 
Fig. ||). Increasing the intensity of the ISRF produces an 
increase in the dust temperature which is felt by the gas 
only in the central parts of the cloud, where the coupling 
of dust and gas is more efficient. As a result, the cen- 
tral gas temperature T g c increases slightly, whereas the 
gas temperature at the boundary T g .b remains unchanged 
(Tg.b — 10 K), since for densities n ~ 10 3 cm" 3 both the 
gas cooling rate and the the cosmic-ray heating rate are 
roughly proportional to n. Thus, one can extend Eq. (^) 
to polytropic clouds bounded by the same external pres- 
sure and characterized by the same thermal sound speed 
ab at the boundary, and write 



Mr, 



Pext 

Table [| shows the value of a n and the density contrast 
Pc/ Pb for clouds with different polytropic index n. For de- 
creasing positive n, the coefficient a n increases, indicating 
a larger critical mass than an isothermal sphere bounded 
by the same external pressure and with the same tempera- 
ture at the boundary. For n — 6, for example, the critical 
mass is expected to be a factor 2.06/1.18 ~ 1.75 larger 
than critical Bonnor-Ebert mass for the same gas temper- 
ature at the boundary (Mbe — 2.6 M for T g = 10 K, see 
Eq.|). 

7. Discussion and Conclusions 

We have examined the gas temperature distribution to 
be expected in interstellar pre-protostellar cores heated 
by the external ISRF. We find that when (as in observed 
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mainly because cooling by gas-grain collisions becomes 
dominant. 

The observed values of the temperature of around 10 K 
in many pre-stellar cores allow limits to be placed on the 
cosmic ray ionization rate similar to the standard value of 
order 10" 17 s" 1 based on the measured cosmic ray flux. 
The fact that measured gas temperatures in cores show so 



little spatial variation (Benson & Myers 1989 



Tafalla et 



central H 2 density n c (cm 3 ) 



Fig. 6. Instability in model cloud cores for isothermal 
and non-isothermal models. The cloud's mass is plot- 
ted as i function of the central H2 donoity for p (iitt — 



2 x 10 " K cm' 3 {solid curves, otablo equilibria; daahed 
curves, unstable equilibria) . Instability sets in at the max- 
imum of M(n c ). The curves are labelled with the values 
of Go- The dotted lines separate the regions of stability 
and instability for this value of p ex t and arbitrary values 
of the temperature (in the isothermal case) and Go (in the 
non- isothermal case) . The thick curve show the mass ver- 
sus central density relation for an isothermal sphere with 
T g = 10 K. 

Table 3. Critical mass and density contrast as 

FUNCTION OF POLYTROPIC INDEX 
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20 


1.39 


15 


00 
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cores), the central density exceeds 3 x 10 4 cm" 3 , there 
is coupling between the gas and dust temperatures and 
hence the gas temperature (like the dust) decreases with 
decreasing radius. At larger radii and smaller densities, the 
dust and gas decouple and the gas temperature may (for 
low external pressures as in the Taurus cloud) decrease to- 
wards the values expected for heating by galactic cosmic 
rays of around 10 K. The region where gas and dust tem- 
peratures are coupled is somewhat interior to the region 
where CO is highly depleted. However, CO depletion does 
not seem greatly to affect the temperature distribution 



al. 2002) suggests to us that the observed thermometers 
mainly trace layers of moderate depletion where gas-grain 
coupling is not playing an important role. 

We have also examined the consequences of such a 
temperature distribution for the density distribution in 
hydrostatic equilibrium cores. The changes caused by the 
"real temperature distribution" are minor and the char- 
acteristics of a "marginally stable Bonnor-Ebert sphere" 
are similar to those in the isothermal case. 

An interesting point which emerges from these calcu- 
lations is that the temperature in the core nucleus is sen- 
sitive to the external radiation field (because the core nu- 
cleus as a rule is heated by small particle MIR emission 
from the borders of the surrounding cloud). This is in par- 
ticular the case for cores with high external pressure (and 
hence high density) such as those studied by Johnstone 
ct al. (200l| , |2000D in the Orion B and pOph clouds where 
Pext appears to be above 10 6 cm _3 K. In these cases, we 
expect the dust temperature to follow the gas tempera- 
ture and thus the behavior should be roughly like a neg- 
ative index polytrope. The external radiation fields are 
also higher in Orion and Ophiuchus than in Taurus and 
thus the core masses may also rise somewhat. For con- 
stant external pressure, the mass of the marginally stable 
Bonnor-Ebert sphere rises with about the 0.35 power of 
the external radiation field. It would be useful to have di- 
rect temperature estimates to confirm these expectations. 

Another result of this study is that high mass ther- 
mally supported cores are incompatible with high col- 
umn density. Thus it is difficult to imagine observed 
high column density (greater than 10 22 cm' 2 ) high mass 
(greater than 10 M Q ) cores going through a series of quasi- 
equilibrium states en route to collapse. Hence we suspect 
that higher mass cores are either not stable or have other 
(magnetic) means of support. These other means of sup- 
port become apparent observationally both because ob- 
served cores show large departures from spherical sym- 
metry and also because the observed line widths in most 
cores can only be explained in terms of turbulence and 
often of supersonic turbulence. 

Thus thermally supported cores may be the exception 
rather than the rule. It is significant nonetheless that cores 
in regions such as Taurus where predominantly low mass 
star formation is taking place have in general line profiles 
showing a large component of thermal broadening. This 
suggests that where (as in Taurus), the pivotal state may 
be a core with predominantly thermal support, only low 
mass stars are likely to form. In clouds such as Ophiuchus 
and Orion with higher radiation fields and pressures and 
with cores having predominantly non-thermal support, 
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higher mass stars may become possible and star forma- 
tion may take another course. 

It is also significant that occasionally, one finds cases 
like B68 where the data are consistent with hydrostatic 
equilibrium (marginally stable) and pure thermal support. 
In fact, we find that when one takes the gas temperature 
dependence into account, B68 is marginally stable. But 
irrespective of whether this is true or not, B68 gives every 
sign of being close to the pivotal state from which proto- 
stellar collapse will commence. 
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